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Quantum Corrections in Quintessence Models 
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We investigate the impact of quantum fluctuations on a liglit rolling quintessence field from three 
different sources, namely, from a coupling to the standard model and dark matter, from its self- 
couplings and from its coupling to gravity. We derive bounds for time-varying masses from the 
change of vacuum energy, finding Ame/rue <C 10~^^ for the electron and Anip/mp <C 10~^^ for the 
proton since redshift 2 ~ 2, whereas the neutrino masses could change of order one. Mass-varying 
dark matter is also constrained. Next, the self-interactions are investigated. For inverse power 
law potentials, the effective potential does not become infinitely large at small field values, but 
saturates at a finite maximal value. We discuss implications for cosmology. Finally, we show that 
one-loop corrections induce non-minimal gravitational couplings involving arbitrarily high powers 
of the curvature scalar R, indicating that quintessence entails modified gravity effects. 

PACS numbers: 95.36.+X, 04.62.-|-v, ll.lO.Gh 



I. INTRODUCTION 

A possible explanation for the observed [1 , 2] accelera- 
tion of the universe is given by a light rolling scalar field 
[3, 4], usually called quintessence field. The dynamics of 
this field can lead to a decaying dark energy, and thus 
address the question why the cosmological "constant" is 
small but non-zero today. Presently, there is a huge num- 
ber of dark energy models based on scalar fields, see e.g. 
ref. [5] for a review. Models admitting e.g. tracking 
solutions [6] or general scaling solutions [5] additionally 
possess some appealing properties like attractors which 
wipe out the dependence on the initial conditions of the 
field in the early universe, or a dynamical mechanism nat- 
urally yielding an extremely small classical mass of the 
quintessence field of the order of the Hubble parameter. 
The latter can be necessary e.g. to inhibit the growth of 
inhomogeneities of the scalar field [4] . 

However, the rolling quintessence field usually cannot 
be regarded as completely independent of other degrees 
of freedom. If the quintessence dynamics are for exam- 
ple governed by a low-energy effective theory which is 
determined by integrating out some unknown high en- 
ergy degrees of freedom, involving e.g. quantum gravity, 
string theory or supergravity [5, 7], the low-energy theory 
should generically contain couplings and self-couplings of 
the quintessence field suppressed by some large scale, e.g. 
the Planck scale. In some cases such couplings can be di- 
rectly constrained observationally, like for a coupling to 
standard model gauge fields [8] , whereas a significant in- 
teraction with dark matter seems to be possible [9] and 
is used in many models, e.g. [10-15], often accompanied 
by a varying dark matter mass (VAMP) [16-19]. Models 
leading to time- varying standard model masses and cou- 
plings, including mass- varying neutrinos (MaVaNs), from 
a corresponding coupling to a rolling dark energy field 
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are also frequently considered, see e.g. [20-29], which 
can lead to potentially observable effects like a variation 
of the electron to proton mass ratio [30, 31] and the fine- 
structure constant [32, 33] or violations of the equiva- 
lence principle [25, 34], and could have an effect on BBN 
[34, 35]. Furthermore, higher order self-interactions of 
the field seem to be a typical feature necessary for a suc- 
cessful dark energy model, involving e.g. exponentials or 
inverse powers of the field [3, 4, 6, 36-38]. Non-minimal 
gravitational couplings of the scalar field have also been 
studied in various settings [39-46], constrained e.g. by 
solar system tests of gravity and BBN. 

Because of the presence of quantum fiuctuations of the 
standard model degrees of freedom as well as of dark mat- 
ter and of dark energy itself, the dynamics of the quint- 
essence field receive radiative corrections. Therefore, it 
is important to study the robustness against these cor- 
rections, see e.g. refs. [22, 47-52] for previous work. 
Apart from the long-standing problem of the overall nor- 
malization of the effective quintessence potential (i. e. 
the "cosmological constant problem"), which is not ad- 
dressed here, quantum corrections can influence the dy- 
namics e.g. by distorting the shape or the flatness of the 
scalar potential. In the case of a coupling to standard 
model or heavy cold dark matter particles, this leads to 
tight upper bounds for the corresponding couplings, due 
to a physically relevant field-dependent shift in the corre- 
sponding contribution to the vacuum energy [53, 54]. The 
quantitative bounds obtained in this way can be trans- 
lated into bounds for time-varying masses and for the 
coupling strength to a long-range fifth force mediated by 
the quintessence field, as will be shown in section II. The 
self-coupling and gravitational coupling of the dark en- 
ergy field are necessary ingredients of basically any given 
model. The corresponding quantum corrections can lead 
to significant modifications of the shape of the scalar po- 
tential, discussed in section HI using an infinite resum- 
mation of bubble diagrams. This accounts for the fact 
that the effective theory for a scalar field does not de- 
couple from the high energy regime due to the presence 
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of quadratic divergences. We also discuss implications 
for cosmology. In section IV we investigate which kind 
of non-minimal gravitational couplings are induced by 
quantum fluctuations of the dark energy scalar field. 



II. LIFTING OF THE SCALAR POTENTIAL 

Generically, the light mass of the quintessence field 
is unprotected against huge corrections induced by the 
quantum fiuctuations of heavier degrees of freedom. Fur- 
thermore, not only the mass but also the total potential 
energy have to be kept small, which is the "old cosmo- 
logical constant problem". In General, the explanation 
of the present acceleration and the question about huge 
quantum field theoretic contributions to the cosmological 
constant may have independent solutions. However, it is 
required in quintessence models that the total cosmolog- 
ical constant is small. 

Here we will take the following attitude: Even if we ac- 
cept a huge amount of fine-tuning and choose the quint- 
essence potential energy and mass to have the required 
values today by a suitable renormalization, there may be 
huge corrections to the potential at a value of the quintes- 
sence field which is slightly displaced from todays value. 
Since the scalar field is rolling, such corrections would 
aflFect the behaviour of the quintessence field in the past, 
and could destroy some of the desired features (like track- 
ing behaviour) of dynamical dark energy. 

To calculate the effect of quantum fluctuations, we will 
impose suitable renormalization conditions for the effec- 
tive quintessence potential. Therefore, we are going to 
argue that under certain general prerequisitions, there 
remain only three free parameters (linked to the quar- 
tic, quadratic and logarithmic divergences) that can be 
used to flx (or fine-tune) the effective potential which is 
induced by the fluctuations of heavier particles coupled 
to the quintessence held. Following the above argumen- 
tation, these free parameters will be fixed at one-loop 
level by imposing the renormalization condition that the 
quantum contributions to the effective potential and its 
first and second derivative vanish today (0o = (t^ito))' 

V{cP^<Po)iL = 

V'{ct> = <l>o)iL = (1) 

V"{(^ = <f>0)lL = 0, 

where V{(I))il denotes the one-loop contribution to the 
effective potential Vejj{(j)) = V{<j)) + Vil ((/)) + ... . 

Since the quintessence field generically changes only 
slowly on cosmological time-scales, one expects that the 
leading effect of quantum fluctuations is suppressed by a 
factor of the order 

V"'{<p = <f>o)iL{<p{to)Atf (2) 

(with At of the order of a Hubble time) compared to the 
potential V{(j)). 



The coupling between quintessence and any massive 
particle species j is modeled by assuming a general de- 
pendence of the mass on the quintessence fleld. This 
general form includes many interesting and potentially 
observable possibilities, like a time- varying (electron- or 
proton-) mass mj{<j){t)), a Yukawa coupHng drrij/dcp to 
fermions (e.g. protons and neutrons) mediating a new 
long-range flfth force, or a coupling between dark energy 
and dark matter (dm) of the form (see e.g. [9]) 

d 111 ni,iiii {(t>) ■ 

Pdm + oHpdm = Pdm 0- (3) 

In terms of particle physics, a dependence of the mass 
on the dark energy fleld (j) could be produced in many 
ways, which we just want to mention here. One possibil- 
ity would be a direct (/(-dependence of the Higgs Yukawa 
couplings or of the Higgs VEV. For Majorana neutri- 
nos, the Majorana mass of the right-handed neutrinos 
could depend on (j) leading to varying neutrino masses 
via the seesaw mechanism [55] . The mass of the proton 
and neutron could also vary through a variation of the 
QCD scale, for example induced by a 0-dependencc of 
the GUT scale [56] . Additionally, a variation of the weak 
and electromagnetic gauge couplings could directly lead 
to a variation of the radiative corrections to the masses 
[53] . Possible parameterizations of the dependence are 
to((/)) = mo(l + [3 f {(j) / Mpi)) with a dimensionless cou- 
pling parameter (3 and a function f{x) of order unity or 
m(?!>) =moexp(/J0/Mp;) [22]. 

Induced Effective Potential 

The one-loop contribution to the effective potential for 
the quintessence fleld can be calculated in the standard 
way from the functional determinants of the propagators 
with mass m{(j)): 

) \b (4) 

where B and F run over all bosons and fermions with in- 
ternal degrees of freedom qb and gp respectively, and the 
momentum has been Wick-rotated to euclidean space. To 
implement the renormalization conditions (1), we con- 
sider the class of integrals 

Io{m')^ /(04ln(g^+m^) 

2^= fd'q 1 _ (-l)''-^ 2^ 

^fcl"^ J-y^27r)4(«'+m2)'= " {k-l)\ ((im2)'='°^"' ^' 

(5) 

which are finite for fc > 3. Following the standard pro- 
cedure described e.g. in [57] the divergences in /q, Ii 



3 



and I2 are isolated by integrating 73 with respect to rn?, 
yielding 



/m /""^a /'"^2 



(6) 



+ Dq + Dxrri^ + £>2m^ 



with infinite integration constants Dq, Di and I?2- Thus 
one is led to introduce three counterterms proportional 
to m°, m? and to cancel the divergences, which can 
be easily reabsorbed by a shift of the scalar potential V. 
This leaves a finite part of the same form as (6) but 

with the three infinite constants replaced by three finite 
parameters that have to be fixed by the three renormal- 
ization conditions (1). It is easy to see that the appropri- 
ate choice can be expressed by choosing the lower limits 
in the integration over the mass to be equal to its 
todays value mg: 



/^"'*''(TO2;mg) =2 / dm'i I dm'i I dmi h{mi) 
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(7) 

where 73(771^) = l/{32Tr'^m'^) has been used. 

Thus the renormalized one-loop contribution to the ef- 
fective potential which fulfills the renormalization condi- 
tions (1) is uniquely determined to be 
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(8) 



The effective potential renormalized in this way can be 
regarded as the result of a fine-tuning of the contribu- 
tions from the quantum fluctuations of heavy degrees of 
freedom to the quintessence potential energy, slope and 
mass at its todays values, i.e. evaluated for = ^q. How- 
ever, when the quintessence VEV had different values in 
the cosmic history, the cancellation does not occur any 
more and one expects the huge corrections of order 
to show up again, unless the coupling is extremely weak. 
Indeed, this argument yields extremely strong bounds for 
the variation of the masses with the rolling field (p. Sim- 
ilar considerations have been done e.g. in [53, 54]. To 
obtain a quantitative limit we require that the one-loop 
contribution to the potential should be subdominant dur- 
ing the relevant phases of cosmic history up to now, which 
we take to be of the order of a Hubble time, in order to 
ensure that the quintessence dynamics, e.g. tracking be- 
haviour, are not affected. For the corresponding (/(-values 
this means that we require 

Vil{^) « (9) 

If we Taylor-expand the one-loop effective potential (8) 



around todays value (pQ, the first non- vanishing contribu- 
tion is by construction of third order, 
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dlnm| V"{(\))^ 



(10) 

Here the index j runs over bosons B and Fermions F 

(with the minus sign in front of gj for the latter) , and eq. 
(7) has been used. In the last line, we have rewritten the 
dependence on the quintessence field in a dependence 
on its mass m'^{(j)) = V"{(j)). Typically, the mass is of 
the order of the Hubble parameter, which is today Hq ~ 
lO^^^eV. In many generic scenarios, e.g. for tracking 
quintessence models [6] , the quintessence mass also scales 
proportional to the Hubble parameter H during cosmic 
evolution. Therefore, we assume that 

lnV"{cj>)/V"{<f>o) ~ In 777772 < ^^^^^ ^ ^^^^ 

where z is the cosmic redshift. If we want the inequality 

(9) to hold up to a redshift Zmax^ the most conservative 
assumption is to replace the logarithm in the last line in 

(10) by its maximal value of order 3 ln(l -|- Zmax) and the 
right hand side of (9) by the minimal value V{4>a). Fur- 
thermore, the inequality (9) is certainly fulfilled if each 
individual contribution to the one-loop potential (8) re- 
spects it. Altogether, under these assumptions the re- 
quirement (9) that the quintessence dynamics are un- 
altered up to a redshift z^ax yields the bound for the 
variation of the mass rrij of a species j (with gj inter- 
nal degrees of freedom) with the quintessence mass scale 

V" ~ 772 



dhumj 



dlnV" 



1 



3ln{l + Zmax) ygjinjicpo)"^ 



(12) 



This bound is the main result of this section. It scales 
with mass like m"^/"^, i.e. the bound gets tighter for 
heavier particles. Inserting z^ax ^ ^eq ^ 10^ and ex- 
pressing the potential energy V{(j>o) = ^—j^^de^Q 
terms of the dark energy fraction Q^e and equation of 
state iVde with 77o ^ 70km/sMpc yields 



(13) 



Finally, we want to remark that there remains the pos- 
sibility that several masses rrij change in such a way 
that the total contribution to the effective potential stays 
small [53]. It would be interesting to look for a spe- 
cial (dynamical) mechanism or a symmetry which leads 
to such fine-tuned correlated changes. Otherwise, there 
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seems to be no motivation for such a behaviour. An ex- 
ample for such a mechanism could be based on supersym- 
metry, where the masses of fermions and their superpart- 
ners would have to change in the same way if SUSY was 
unbroken, so that their contributions in eq. (4) would 
always cancel. However, this is not the case below the 
SUSY breaking scale. 



Bounds on Quintessence Couplings 

The upper bound (12) can be directly related to upper 
bounds e.g. for the coupling strength to a long range 
force mediated by the light scalar field, and for cosmic 
mass variation. The relative change of the mass rrij since 
redshift z can be related to the derivative d In rnj /d In V" 
using eq. (11), 



Ar 



rrii 



imj V"((l)) dlnmj 
<31n(l + z)- ' 



d\nV" V"i(j)o) 



d\nV"' 



(14) 



which means the bound (13) directly gives an upper limit 
for the relative mass variation of species j since redshift 
z. For example, for the variation of the electron mass 
since z ^ 2 we find 



Ar 



< 0.7- 10" 



1 — UJde ^de 
2 OJ 



(15) 



which is at least six orders of magnitude below obser- 
vational constraints for a change in the electron-proton 
mass ratio [34]. For heavier particles, the bounds are even 
stronger by a factor (me/m)'*/'^, see figure 1, e.g. of the 
order Amp/nip <C 10~^^ for the proton. The only known 



particles which could have a sizeable mass variation due 
to the bound (13) are neutrinos. Thus models consider- 
ing mass- varying neutrinos and/or a connection between 
dark energy and neutrinos (see e.g. [23, 58]) do not seem 
to be disfavored when considering quantum fluctuations. 
If the bound (13) is saturated, it is even possible that 
backreaction effects could influence the quintessence dy- 
namics in the recent past, where the turnover to a dark 
energy dominated cosmos occurs. 

If we consider fermion masses which depend on the 
quintessence field, the corresponding coupling mediates 
a Yukawa- like interaction with typical range ~ y^i/^ 
Hq^ and Yukawa coupling strength yj = dmj{<f>)/d(f> [4] 
of species j to this fifth force. Inside the horizon, this 
is a long-range interaction like gravity, which could be 
detected via a violation of the equivalence principle. For 
nucleons, this puts strong bounds on the coupling yp^n of 
order 10"^'' [4]. On the other hand, the coupling strength 
is constrained by the bound (12) via the relation 



_ drrij 



1 y'dlnm^ rn 

fYl ■ i_ : 

2 ^ V" dlnV" 2Md\nV"' 



dlnmf 
' ^ (16) 



where we introduced the scale height M = 
{dh\V" / dcj))^"^ of the quintessence mass which is 
typically of the order of the Planck scale [6]. If we 
consider the proton and neutron as effectively massive 
degrees of freedom, we obtain an upper limit from the 
requirement (13) 



V M 



0.7 



(17) 

which is more than ten orders of magnitude below the 
limit from the tests of the equivalence principle [4] , 



typical range for 
bounds from obs., 
e.g. Aoc/a 




10'' 10° 
m/eV 

Figure 1: Bounds for cosmic mass variation since a ~ 2 from 
the radiative correction to the quintessence potential in de- 
pendence of the mass m. The red (vertical) lines mark the 
masses of some standard model particles. The limits inferred 
from observations e.g. of Aa/a strongly depend on the con- 
sidered particle type and further assumptions, but typically 
lie around 10"" to 10"^ [34]. 



Vp.n < 10 
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(18) 



These limits can be compared to the corresponding grav- 
itational coupling ^ rrij/Mpi, e.g. of the order ^ 10~^^ 
for the nucleons. Thus the bound in eq. (13) also directly 
gives a bound for the relative suppression 



dlnm, 



ruj/Mpi d{(j}/Mpi) 



(19) 



of the coupling strength to the fifth force mediated by the 
quintessence field compared to the gravitational coupling, 
giving roughly (for M ~ Mpi, -t- 1 < 1, firfe ~ 0.7) 



Am, / meV\ 
< 4 ' ' 



4/3 



V J 



4/3 

— ) • (20) 
mi- 



Note that the bound from eq. (17) also holds for other 
species (with a scaling ~ m~^^^ with mass), whose quint- 
essence couplings are in general not constrained by the 
tests of the equivalence principle [4]. This is also true 
for dark matter, if it consists of a new heavy species like 
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e.g. a WIMP, which severely constrains any coupling via 

a </)-depcndcnt mass, 

Vdm = dnidm/d^ « 10-'^*^ (TeV/mrf„)i/3^ (21) 

corresponding to a limit of the order 

Amdm/rridm « 10"''' (TeV/md„)^/3 (22) 

for a mass variation between z ^ 2 and now from cq. 
(20). However, this constraint is not applicable if dark 
matter is for example itself given by a scalar condensate, 
e.g. as in axion models. 

III. SELF-INTERACTION OF THE SCALAR 

If the light scalar field responsible for dark energy has 
itself fluctuations described by quantum field theory, its 
self-interactions will also contribute to the effective po- 
tential. Typical potentials used in the context of quint- 
essence, involving e.g. exponentials [3, 4, 36, 37], contain 
self-couplings with an arbitrary number of legs, which are 
suppressed by a scale M, typically of Planck-size. Such 
couplings could arise e.g. as an effective theory by inte- 
grating out some unknown high-energy degrees of free- 
dom. Usually, quantum fluctuations in the presence of 
such couplings can be treated by an expansion in the in- 
verse of the suppression scale M. However, in the case of 
a scalar field, the high-energy sector does not completely 
decouple due to the well-known quadratically divergent 
contributions. In the context of an effective theory, the 
quadratically divergent diagrams, e.g. the tadpole graph, 
are intrinsically governed by a scale A which is charac- 
teristic for the high-energy scale up to which the effective 
theory is valid. In the simplest case, A can be imagined 
as a cutoff for the momentum cycling in the loop. Both 
high-energy scales M and A could of course be related in 
a way depending on the unknown underlying high-energy 
theory. Since the suppression scale M could be as large 
as the Planck scale, it is even possible that the same is 
true for A. However, since unknown quantum gravity ef- 
fects will play an important role in this regime, we will 
just assume an upper bound A < Mpi. 

In order to establish a meaningful approximation, it 
would be desirable to resum all contributions propor- 
tional to powers of A/M, whereas the tiny mass of the 
quintessence field given by V" , which is typically of the 
order of the Hubble scale today, could admit a perturba- 
tive expansion e.g. in powers of V" /A'^. In the following 
we will motivate that such an expansion might indeed 
be possible, and calculate the leading contributions ex- 
plicitely. 

A typical feature of quintessence potentials is that the 
self- couplings V^''^ {(f>) with k lines are suppressed like 



(23) 



with M = {d\nV"/d(t>)-^, where for example M ^ Mpi 
for exponential potentials and M (p for inverse power 



law and tracking potentials [6] with cf) ^ Mpi in the 
present epoch. The effective potential can be calculated 
by the sum over all IPI vacuum diagrams with propaga- 
tor G(p) = i/(p2 - V"{(t>)) and vertices -iV^''\(j)) with k 
legs. If one uses the power-counting estimate in eq. (23), 
then a L-loop diagram with Vk vertices with A; > 3 legs 
and E external lines is (in terms of dimension-full quan- 
tities, i.e. disregarding logarithms etc.) proportional to 



n [y"'') 



^2(L-V+1) 



A2 



/ A^ 



L + l 



M 



A-E 



(24) 

which shows that the diagrams with only one vertex are 
the leading contribution in V" /h? <C 1, whereas there is 
not necessarily a strong suppression of contributions with 
high L. This indicates that in the case of quintessence- 
like potentials an appropriate expansion parameter is 
given by the number of vertices V, whereas the loop ex- 
pansion becomes meaningless if A > M. 



Bubble Approximation 

We will now calculate the effective potential in leading 
order in the number of vertices V and show that it can 
be consistently renormalized for quintessence potentials 
obeying eq. (23), up to higher order corrections. The 
graphs with V = 1 are "multi-bubble" graphs, i. e. graphs 
where an arbitrary number of tadpoles is attached to the 
vertex. Thus, diagrammatically, the effective mass, i.e. 
the second derivative of the effective potential^, is given 
in leading order in V by the infinite sum 



Vi 



oo 

E 

L=0 



F(2i+2)((/.) 



r)4 q2 + v"{cl>) 



m^=V"{(/>) 



(25) 

with euclidean momenta and Ji as defined in (5). The last 
line is a compact notation where the derivatives in the 
exponential act on the tree level mass V"((f)) on the right 
hand side. Note that this is the infinite resummation re- 
ferred to above, in leading order in V. The physically 
relevant scale of the quadratically divergent tadpole in- 
tegral Ii is given by a high energy scale characteristic for 
the effective theory, as discussed above. To implement 
this behaviour in a consistent way, we split 



r — /-finite i rd 
ll=Il +li 



div 



(26) 



^ It is also possible but less convenient to calculate Vejj- directly. 
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where the finite value should represent the physically rel- 
evant part. The formally divergent part can be consis- 
tently reabsorbed by shifting the potential by an analytic 
function of the field value 6, 



V{ct>)^V{ct>) + 5V{ct>), 



(27) 



as is shown in detail in appendix A, up to corrections 
which are suppressed by a factor of the order yy/finitc 
relative to the leading contributions^ and assuming eq. 
(23) for the potential. Thus the bubble approximation 
allows a self-consistent renormalization in leading order 
in y //finite^ with the important result (see appendix A) 



K^(0) = exp 



1 

_ /-finite/ 2n_^ 
2^1 ^™ ^#2 



y"(</.) 



=y"(0) 



(28) 

up to terms which are relatively suppressed by p^'y/fimto 
This result means that the i-loop contribution to V^j^(0), 
see eq. (25), which is proportional to (/i)^, leads to 
a leading-order contribution of the form (/fimto)^ after 
renormalization, i.e. "renormalization and raising to a 
power do commute" in this case. 

The physically relevant part /fi"'*<= will be determined 
by an underlying theory, possibly involving quantum 
gravity effects, which produces the quintessence potential 
as effective theory. In the present approach we param- 
eterize this unknown function by a Taylor expansion in 
V", 



167r2 



(29) 



where the dots include linear and higher terms in V" . 
Since V" is typically of the order of the Hubble scale, 
we neglect these contributions with respect to the high- 
energy scale ji < Mpi , which should maximally be of the 
order of the scale A introduced in the first part of this sec- 
tion. Furthermore, one cannot exclude a priori the possi- 
bility that /fi"''<' can be negative or positive, indicated by 
the two signs. This again depends on the embedding into 
the underlying theory, where unknown quantum gravity 
effects could play a major role. Additionally, there are 
examples like the Casimir effect, where it is known that 
the sign of the rcnormalizcd 00-componcnt of the energy- 
momentum tensor can be negative or positive, depending 
e.g. on boundary conditions and geometry, even though 
the unrenormalized contribution is positive definite. 

Altogether, the leading contribution to the effective 
potential is given by 



Veffi^) = exp ±—^fi 



y(</.), (30) 



^ In order to rcnorrnalize also the sub-leading contributions, it 
would be neeessary to include also two-vertex graphs since they 
are of the same order. 



up to corrections suppressed by V" / <C 1. This is 
the main result of this section. It has been obtained by 
integrating eq. (28) twice with respect to </> after using 
eq. (29) and absorbing the constant of integration, which 
corresponds to the quartic divergence, into 5V. Thus the 
potential is, as usual, only determined up to an arbitrary 
additive constant, which we set to zero, corresponding to 
the unresolved "cosmological constant problem". 



Stability of Quintessence Potentials 

The above result (30) gives a simple prescription to 

estimate the stability of a quintessence potential V{(t>) 
under its self-interactions. In the following, we will in- 
vestigate the effect of an operator of the form 



exp ( ± 



(31) 



on some typical potentials often used in dynamical dark 
energy scenarios. One archetype class of potentials are 
given by (combinations of [59, 60]) exponential poten- 
tials [3, 4, 37]. Remarkably, an exponential of the field <p 
is form-invariant under the action of the operator (31). 
Consider e.g. the following finite or infinite sum of expo- 
nentials, 



y(^)=5:F,exp(-A,^). 



(32) 



The only effect of (31) is a simple rescaling of the pref- 
actors Vj according to 



2 ,,2 



Vj ^ y. exp ± 



Am 



(33) 



This extends the result of ref. [22] for the one-loop case, 
which would corresponds to the first term in a Taylor 
expansion of (31). Note that if ~ Mpi the correction 
can be of an important size, and can influence the rela- 
tive strength of the exponentials in (32). The necessary 
condition of validity V"{<j)) -C for the bubble approx- 
imation is typically fulfilled when V <C ^I'^M^i, which 
implies that it is applicable if /i ^ H^ax: where H^ax is 
the maximum value of the Hubble parameter where the 
field (j) plays a role. For example, H^ax could be the in- 
flationary scale Hinf- e. g. around IQ-'^'^GcV. Altogether, 
exponentials seem to be stable under the considered ra- 
diative corrections. 

Another often discussed class of potentials are (combi- 
nations) of inverse powers of the field (j) [4, 6, 22, 48], 



y(.^) = ^c„<i 



(34) 
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The action of the operator (31) yields 

(35) 

where the F-function inside the sum over L has been 
replaced by its definition via an integration over the pos- 
itive real axis in the second line. This integral only gives 
a finite result if the negative sign in the exponent is used, 
which we will therefore assume from now on. We will first 
discuss two limiting cases where the integral can be easily 
solved analytically. For large field values <j> ^ which 
corresponds to small potential energy and -curvature, the 
second term in the exponent — t — /i^/(327r^0^)t^ appear- 
ing in (35) can be neglected, which implies that asymp- 
totically 

v;jy(0)^y(^) = ^c„r", (^^oo. (36) 

a 

This means the low energy regime where V and its deriva- 
tives go to zero is not changed. For the opposite limit 
where cj) <^ fi, the integral in the last line of (35) can be 
calculated by neglecting the first term in the argument 
of the exponential, 

v-r(f) A* V" t 

= > tt^tt-tCq ^ ~ const. 

^ 2F(a) U7TV2J 

Thus the effective potential approaches a constant finite 
value for 4> < yLt/(47r\/2) of the order V{fj,), which gets 
smaller for larger values of fi, see figure 2 for the special 
case V oc (f>^^. Furthermore, it is easy to see that also 
V^y{(j)) approaches a constant value 

K>(«-i:^=.(^) ^ (3S) 

These results show that the singular behaviour of the 
potential V{4i), see eq. (34), for — > is not present 
in the effective potential, where a constant value of the 
order V{fi) is approached instead. 

The bubble approximation requires that V" <C /i^, 
which means we again find that the approximation is 
valid as long as » Hmax, with Hmax < H,nf, as for 
the exponential potential, if we simply assume tracking 
behaviour [6]. If the scale /i is very large, the bubble 
approximation is always valid in the range of cosmolog- 
ical interest anyway. For example, for"^ V oc (/)~^ as in 



For V = Ca(p " typical values within quintessence scenarios are 
~ H^lVr+^, e.g. C2 ~ (100MeV)6 [6]. 



figure 2, one finds that the condition V" <^ /i^ with 
fi > 10~^Mpi is always fulfilled as long as 4>/<j>o > 10^^®, 
where c/jq Mpi is the field value today, which is far 
below the relevant range of (j). Furthermore, eq. (38) 
shows that due to the absence of a singular behaviour in 
the second derivative of the effective potential the condi- 
tion V^'ff/fi^ < 1 can hold for all positive values of (j) if 

/i ^ maxQ,(cQ,)^/^"+'''. This might indicate that the bub- 
ble approximation is indeed applicable in the total range 
of field values. However, to show this formally it is nec- 
essary to additionally perform a 2PI approximation [61] 
where a dressed propagator containing V"jj instead of V" 
appears in the loops. This is left to future work. 




10"^ 10"^ 10"'* 10"^ 10'^ 10"'' 1 



Figure 2: Effective potential Veff{(j)) for an inverse power law 
potential V{(j)) oc (f>~°' with a = 2 for various values of fi 
on a double logarithmic scale and normalized to todays value 
V{(l>o)- From top to bottom, /i is enlarged by a factor 2 for 
each red (thin) line. The black (thick) line is the tree potential 
V{4')- For ^ /i/(47r\/2) the effective potential approaches a 
constant value, whereas V gets huge. The redshift-scale on the 
right-hand side applies for the tracking solution of V only and 
illustrates when the deviations of Veg from V become relevant 
in cosmic history going backward from (j)/4'o ~ 1 (today). 

Let us now estimate in how far typical tracking quint- 
essence models are changed by considering the effective 
potential from eq. (35). Since the field value today is typ- 
ically of the order of the Planck scale [6], the large-field 
limit eq. (36), where the effective potential approaches 
the tree level potential and the corrections are negligi- 
ble, is only applicable when /i <^ Mpi. For values up 
to /i < Mpi/ 10 the field cj) can have a tracking solution. 
The redshift Equant in cosmic history where the effective 
potential starts to deviate from the tracking potential, 
see figure 2, gives a rough estimate at which redshift the 
tracking sets in. For a potential dominated by a single 
inverse power V oc 4>~" we obtain, requiring a deviation 
of the effective potential of less than 1% and using the 
tracking solution during matter and radiation domina- 
tion with equation of state uj^ — ^j^(l + ojb) — 1 [6], 
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with ojb = 0, 1/3 respectively, 



-0.5 



•2^quant 



/i/(4^V2«(a + 1)) 



e.g. assuming fj, ~ Mpi /lOO [Mpi 



Equant 



300 for a = 2 and z, 



quant 



(39) 



1/vG) one gets 
130 for a = 1. 



Similar bounds should also hold for other types of po- 
tentials, e.g. like the SUGRA-potential [48], which are 
dominated by an inverse-power law behaviour at red- 
shifts z ^ 0.5. For values /i > Afp//10, there can be 
large deviations from the tracking solution even at low 
redshifts and today, see figure 3 for an exemplary case 
with V oc (/)"". If jjL is extremely large, there is a direct 
transition from the slow roll regime with 4> fJ., equa- 
tion of state ~ —1 and dark energy fraction fi^ <^ 1 
in the flattened effective potential Veff to the Dark En- 
ergy dominated accelerating solution for (j) > AIpi with 
51^ — > 1 and —1, and thus the solution never per- 

forms scaling with uj^ — — as for V. In the case 
a = 1, the equation of state today tode = uJip{z = 0) is 
enhanced for 0.1 < fJ^/Mpi < 1.3 compared to the track- 
ing value, and gets smaller for even larger^ fi, see figure 
3. Moreover, the sign of du^/dz can change depending 
on ^. 



IV. EFFECTIVE ACTION AND CURVED 
BACKGROUND 

Since any dynamical dark energy scenario is generically 
used in a curved space-time setting, e. g. described by 
a Robertson- Walker metric, it is important to study also 
the quantum fluctuation on such a background. Gener- 
ically, dynamical dark energy scenarios making use of a 
scalar field involve non-renormalizable interactions sup- 
pressed by some high-energy scale up to the Planck scale. 
Therefore it is important to include such interactions 
when discussing quantum induced non-minimal couplings 
between the dark energy scalar field and gravity. 

In this section, we will investigate the behaviour of a 
scalar field whose potential contains at least one term 
which is non-renormalizable in the usual power-counting 
sense, e. g. V{(t)) 9 0" with n > 4, using the semiclassical 
treatment of curved space-time in quasi-local approxima- 
tion. This includes most quintessence models as well as 
models that can be reformulated using a canonical scalar 
field. 

On a Minkowski background, the one-loop contribu- 
tion to the effective potential requires the introduction 
of counterterms proportional to V"{(t>) and If 
V[(j)) is an analytic function of and contains at least 




0.2 0.4 0.6 0.8 




* Note that even when > Mpi the renormahzed tadpole 7^"'*'= 
can still be sub-Planckian due to the loop factor I/IGtt^. 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
a 



Figure 3: Top: Evolution in the (O^, a;0)-plane for the effec- 
tive potential Veff of an inverse power law potential V oc 
with Q = 1 for various values of n keeping Hq = 73km/sMpc 
and ^Ide = i^4,{z = 0) = 0.76 fixed. From top to bottom, fi 
is enlarged by 0.2A4pi for each red (thin) line starting from 
fi = 0.1 Mpi. The black (thick) line is the tracking solution 
in the tree potential V{(j)), from which the solutions deviate 
considerably for large /i. The four arrows on each trajectory 
mark the points with redshifts z — 2,1,0.5,0.1 from left to 
right. 

Bottom: Contour plot of the equation of state ujdc today 
(2 — 0) using the effective potential V^g of a potential 
V oc (j)'" depending on /i and a. fi — corresponds to the 
case Veff = V. Again, we fixed Ho and ^Ide ~ 0.76 as above. 



one power of the field larger than 4, this immediately im- 
plies that in order to remove all divergences it is necessary 
that V{(j)) contains terms with all powers of (j). In other 
words, the structure of the one- loop quantum corrections 
enforces to keep either only renormalizable terms in Vi^cj)) 
or to admit an arbitrary analytic form for V{(j)). 

In the case of a renormalizable potential in curved 
space-time it is well known [62-64] that in order to be 
able to remove all divergences at the one-loop level it 
is necessary to include also non-minimal coupling terms 
with the curvature scalar R of the form (j>'^R as well as 
terms proportional to R, R^ , the square of the Weyl ten- 
sor C = R^,^paR'"''"'-2R^^Rt"' + ^R'^, the Gauss-Bonnet 
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invariant G = Rf,^,paR^"'P'^-ARf,^,R^"'+R^, OR and 
where □ is the covariant D'Alembertian. The latter three 
terms are total derivatives and thus not relevant for the 
dynamics, but they are needed for the cancellation of di- 
vergences and do appear in the dynamics if their running 
is considered [63]. 

The question we want to address here is how this set of 
operators is enlarged when also non-renormalizable terms 
appear. This means we want to find an (infinite) set 
of operators that closes under one-loop quantum correc- 
tions, i.e. any divergence produced by these operators 
can be reabsorbed into one of them. Therefore we use 
a quasi-local approximation accessible through the heat 
kernel expansion in the form suggested by [65] . The out- 
come of this approach is that, in addition to the operators 
C and G needed in the renormalizable case, it is nec- 
essary to use a generalized potential 1^(0, R) that is an 
arbitrary analytic function of the field (j) and of the curva- 
ture scalar R and furthermore to include total derivative 
terms of the form DB{(f), R) where B{(j), R) is also an ar- 
bitrary analytic function of 4> and R. In other words, the 
"smallest" Lagrangian which is stable under one- loop cor- 
rections in curved space and includes non-renormalizable 
interactions is of the form 

^=1 {d(j>f - V{<i}, i?) + eiC + e^G + □S((/., R), (40) 

with dimensionless coupling constants ei and e2, and ad- 
mitting arbitrarily high powers of both <j) and R to be 
contained in V and B. 

We will now show this result. In Appendix B, where 
we also calculate the one-loop effective action explicitely 
using zeta-function regularization, we show that apart 
from counterterms proportional to C and G one has to 
introduce a counterterm proportional to 



(41) 



i.e. all operators contained in this expression have to be 
present in the tree level Lagrangian. Now assume for 
the moment that V is an (arbitrary) analytic function 
/o(0) only of (j), i.e. can be written as a series in 
with n > 0, as required in flat space. The Lagrangian 
then has to include all terms contained in (/q — R/6)^, 
especially those proportional to (/)"~2/j for all n > 2. 
Consequently, one has to add a term of the form fi((j))R 
to the Lagrangian with an (arbitrary) analytic function 
/i ((/)), giving the counterterm (41) the form (/g + fi'R — 
i?/6)2. This implies the presence of R"^ -terms for all 
n > 2 in the Lagrangian, i.e. a term f2{<j))R^ also has to 
be added, and so on. Recursively, this implies that V = 
X^n fn{<l>)R^, i-e. V is an arbitrary analytic function of 4> 
and R, as was to be shown. Furthermore, counterterms 
proportional to 0(d'^V{(t), R) / dcf)^ — R/5) (see appendix 
B) immediately imply that also B = B{(p, R). 

Some comments arc in order: 
(i) The upper result can also be rephrased by stating that 



even if one starts with a tree level Lagrangian where the 
scalar is minimally coupled to gravity at a certain scale, 
quantum fluctuations will induce non-minimal coupling 
terms with arbitrarily high powers of R in the effective ac- 
tion through the running of the corresponding couplings, 
(ii) There arc two special cases where the introduction 
of arbitrarily high powers of R can be avoided in a way 
which is stable under one-loop corrections: First, as ex- 
pected, the renormalizable case, where for the choice 
V{<l), R)= A + m202/2 + -t- ^(l)'^R/2 + M'^R + cqR'^ 
and B(0, R) = esi? -|- 77^2 all divergences can be ab- 
sorbed, as already mentioned above. Second, the con- 
formally coupled case, where V{(f),R) = V{(j)) +^(j)'^R/2 
with ^ = 1/6, and B{(p,R) = eai? + B{(j)). In this case, 
all divergences proportional to higher powers of R van- 
ish, since R is canceled in eq. (41). In other words, 



if we write V 



y y 



^R"^, there is a fixed 



point where all couplings Cnm for m > 1 do not run for 
C21 = C/2 = 1/12, Cnm = for n > 2, rn > 1 and arbi- 
trary'' cim and com for m > 1. However, within many 
typical quintessence scenarios where the field value 4> is 
of the order of the Planck scale today [6] , a non-minimal 
coupling of the form S^(jpR/2 is restricted since it leads 
e.g. effectively to a variation of the Newton constant 
(167rGe/j)"^ = (IOttG)-! - ^^72. Limits foimd by sev- 
eral authors [39, 40] for specific models lie in the range 
1^1 < 3 • 10~2 which is far below the conformal coupling, 
(iii) It is possible to rewrite the Lagrangian (40) so that 
R enters only linearly by performing a conformal trans- 
formation 9,j,y e^'^g^iu to the Einstein frame. However, 
this introduces an additional scalar a coupled to (j) into 
the theory. Technically, one rewrites the Lagrangian (40) 
in the form 



C = -{d^f-V{<j>,A)-B{R-A)+£„ (42) 

with auxiliary fields A and B, see e.g. [66]. Cc denotes all 
additional contributions including the terms proportional 
C and G. After eliminating B = dV/ dA and performing 
the conformal transformation with a = ln{2K'^dV/dA) 
(k;2 = SttG) the corresponding Lagrangian is 



^-2l + 2^(^^)^ 



(43) 

Apart from the cr-dependent kinetic term of (f> both 
fields interact through the potential given hy U = {V — 
AdV/dA)e~'^'^, where A = A{a,(j)) is given by the in- 
version of (7 = \n(2K'^dV/dA) w.r.t. A. However, the 
physical equivalence of conformally related frames is not 
manifestly obvious (see discussion in [66]) and the calcu- 
lation of quantum corrections in the different frames can 
yield inequivalent results due to a nontrivial Jacobian in 



Since these couplings do not receive any quantum corrections in 
the conformally coupled case, they could be set to zero by hand. 
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the path integral [22]. Thus the main clue is that the non- 
minimal coupling with higher powers of R of the general 
form V{(j), R) cannot be simply rescaled away without 
profoundly changing the scalar sector of the Lagrangian. 
(iv) Since the potential is non-renormalizable it is neces- 
sary to introduce the infinite set (^"i?™ (with n, m > 0) of 
operators to cancel the one-loop divergences. This means 
that the corresponding couplings are not predicted by the 
theory, at least at a certain reference scale, but have to be 
determined in principle by comparison with experiment. 
Of course, this is far from being possible. Nevertheless, 
the result suggests that the framework for searching for 
an explanation of cosmic acceleration could be a combi- 
nation of the two extreme cases of quintessence models 
with V = V{(f>) on the one hand and modified gravity 
scenarios corresponding to V 3 f{R) (see e.g. [5, 67, 68] 
for reviews) on the other. 

The main result of this section is that whenever one 
considers a scalar with non-renormalizable interactions 
and non-conformal coupling, one-loop quantum fluctua- 
tions will induce the presence of terms with arbitrarily 
high powers of the curvature scalar R in the action. In 
the case of dark energy, this means that the quantum 
fluctuations in the quintessence field could lead to a modi- 
fied gravity which differs considerably from the standard 
Friedmannian behaviour. Consequently, in a qiiantized 
picture a scalar condensate with non-renormalizable po- 
tential accounting for dark energy goes hand in hand with 
a modified gravity theory described by a generalized po- 
tential V{(j),R). In fact, such a potential could not only 
give rise to an explanation of the accelerated expansion 
due to a quintessence-like behaviour, but also through 
the modification of gravity. 



V. SUMMARY AND CONCLUSIONS 

In this work we have investigated the efi^ect of quan- 
tum fluctuations in the context of typical dynamical dark 
energy scenarios like quintessence models, where an ex- 
tremely light rolling scalar field supplies the present cos- 
mic acceleration, in some sense similar to the inflaton in 
the early universe. 

First couplings between the quintessence fleld and 
heavier degrees of freedom, like the standard model 
fermions or dark matter, have been discussed. We con- 
strained the discussion to couplings that can effectively 
be written as a field-dependent mass term. These cou- 
plings have to be extremely small even though we fine- 
tune the energy density, slope and mass of the quintes- 
sence field at its todays value by appropriate renormaliza- 
tion conditions for the quartic, quadratic and logarithmic 
divergences in the induced eff'ectivc potential. This leads 
to a bound on time-varying masses between z ^ 2 and 
now of the order 10~^^ for the electron and scaling pro- 
portional m~^/^ with mass, assuming the mass variations 
are not themselves finely tuned in such a way that the 
total shift in vacuum energy is negligible. Moreover, we 



found that the coupling strength to a fifth force mediated 
by the quintessence field has to be suppressed by a num- 
ber of the same order relative to its gravitational coupling 
strength. Only neutrinos could have a large mass varia- 
tion and interact with the quintessence field as strong as 
with gravity. 

Second we introduced a suitable approximation scheme 
to investigate the impact of quintessence self-couplings 
on the shape of the effective potential, while an unde- 
termined additive constant has been fine-tuned to be 
zero, thus bypassing the unresolved "cosmological con- 
stant problem". We showed that the quantum fluctua- 
tions to the scalar potential can be consistently renor- 
malized in leading order in V"/is^, where is a high 
energy scale characteristic for an underlying theory and 
V" the square of the quintessence mass assumed to be 
of the order of the Hubble parameter. While potentials 
involving exponentials just get rescaled, inverse power 
law potentials are flattened at small field values. The ef- 
fective potential approaches a finite maximal value, thus 
truncating the singular behaviour of the inverse power 
law in the field range of interest. This distortion of the 
potential can directly play a role cosmologically if /i is 
large, roughly /i > Mp// 10, and moreover affect general 
properties like tracking behaviour. 

Third we investigated non-minimal gravitational cou- 
plings induced by quantum fiuctuations. Since quintes- 
sence potentials usually cannot be taken to be a quar- 
tic polynomial in the field, non-minimal couplings apart 
from a term of the form (jx^R can be induced. We showed 
that at one-loop all couplings of the form (f>"R'^ with in- 
tegers n and m have to be included, and will be induced 
by quantum corrections unless the field is exactly con- 
formally coupled. Moreover, we showed that this type 
of non-minimal coupling which is nonlinear in R cannot 
be simply removed by a Weyl rescaling, but corresponds 
to a theory with two interacting scalars in the Einstein 
frame. Altogether, this may indicate that the origin of 
cosmic acceleration could not purely be the effect of one 
rolling scalar field, but also involves modified gravity ef- 
fects similar to /(_R)-theories when quantum fluctuations 
of the scalar are taken into account. 

We conclude that quantum fluctuations do play an im- 
portant role (i) in coupled dynamical dark energy scenar- 
ios even if we allow fine-tuning in the form of rcnormal- 
ization conditions, (ii) for the shape of the quintessence 
potential and (iii) for its interplay with gravity. 



Appendix A: RENORMALIZATION IN 
BUBBLE- APPROXIMATION 

In section III the effective scalar potential has been 
calculated using a "multi-bubble" approximation where 

only diagrams with one vertex but with arbitrary num- 
ber of loops have been kept. Here we will show that 
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this effective potential can be consistently renormalizcd 
under the following assumptions: (i) The higher deriva- 
tives of the potential V{(j)) can be roughly estimated 
by = 0{V"{(f>)/M{(f>)'') with the scale height 

M{(p) = {d\nV"/d(f)y^. (ii) The renormahzed value of 
the one- loop tadpole integral /i, see cq. (5), can be 
characterized by a scale n, i. e. Jj^"!*" _ 0(^2-j^ with 
» V"{(f>). (iii) The renormahzation is carried out 
only up to contributions proportional to V"{(f))/ <C 1. 

The third assumption is necessary since the corrections 
to the bubble approximation from graphs with at least 
two vertices are generically also of the order V" {(f) / ji^ , 
and thus would have to be included if terms of this order 
were to be considered. 

The counterterms will in general contain the divergent 
parts of the one- loop integrals Ik- The splitting 



— ^i. 



ir for fc = 0.1,2 



into finite and divergent parts will depend on the reg- 
ularization which is used and a specific renormahza- 
tion scheme. Here, we will not pick up a special pre- 
scription, but perform the renormahzation in full gen- 
erality except the assumption that both the finite and 
divergent parts separately obey the relation //c(rn^) = 
(-1)'=-V(fc - l)\{d/dm^)''Io{m'^), see eq. (5), where it 
is understood that /^'^ = for fc > 3. This implies that 
/^'^(m^) for k = 0, 1,2 are polynomials in m? = V" of 
order fc, exactly as required (see also sec. II). In general, 
the divergent parts will depend on some rcgularization 
parameter and go to infinity if the rcgularization is re- 
moved. For example, in the case of a cutoff A one has 
If'' ~ A^, ~ A2 and 7^^ ~ InA^. The corrections 
to the multi-bubble approximation then contain terms 
of the order F'7/f'"''= - V" / ^i^ and V"/lf^' - V"/K^, 
where the latter contribution is suppressed if the cutoff 
is sent to infinity. Thus, V" / /z^ is an upper hound to the 
corrections which we will therefore use throughout the 
calculations for simplicity, formally corresponding to the 
power-counting rule^ 7i ~ 0(/x^). 

For the renormahzation procedure we start with a 
canonically normalized scalar field with bare potential 
Vb(</>) = V{(p) + 5V{(j)), which is split up into the renor- 
mahzed potential V{(j)) and the counterterms 5V{(j)). 
Since there appears no anomalous dimension in the multi- 
bubble approximation, wc do not introduce a correspond- 
ing (unnecessary) counterterm for simplicity. Further- 
more, we split the counterterms into a series 



(A2) 



in powers of an order parameter ?i, which will in the end 
be set to one, together with the replacement d^q — > h d^q 



This rule can of course only be used to determine whether cor- 
rections suppressed by some power of I\ are sub-leading or not, 
but not in the actual leading-order calculation. 



in the loop integrals. Thus we obtain from eq. (25) 



yeff{<t>) = X, \Jj^ 



n 



N 



)4 q2_^_Y>,j^§yf> 

(A3) 

The main task is to expand this expression in H using eq. 
(A2). Diagrammatically, with lines corresponding to an 
(euclidean) propagator l/iq"^ + V"{(j)j), V^^ is equal to 
the bubble sum as in eq. (25) where also diagrams with 
counterterm-vertices arc added and all propagators can 
carry an arbitrary number of (5y'-insertions, denoted by 
crosses with two or more legs respectively. 

The renormahzation can now be carried out order by 
order in h. The contribution of order h} ("one-loop") is 



(Al) v:ff{<p). 



.a 



^ = 



y(4) 



h{V") + dVl'. (A4) 



This fully determines the one-loop counterterm 



chosen such that 



^ — -'i 



y(4)(^) 



-finite/T/// 



(A5) 



(A6) 



At two-loop order, i.e. at order ti^, one gets 




(A7) 

where h = Ik{V"). The counterterm ^v/"^ = £tSV(' 
can be calculated from eq. (A5), 



(4) 



7-div 

rdiv 



1/(6) 



rdiv 



l + O 




where we have used the assumptions and power-counting 
rules discussed above. The main contribution comes from 
the part where d'^/dcj)'^ only acts on the coupling !/('*) in 
eq. (A5), whereas the other terms are suppressed by a 
relative factor of the order V" / fi^ <C 1. The logarith- 
mic divergence 72 belongs to the sub-leading order and 
can only be consistently renormahzed together with two- 
vertex graphs as discussed above. For brevity, we will 
denote any such corrections by e = V"p{l2^'^)/ jj? where 
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p stands for a polynomial with maximally order one co- 
efficients. 

The third diagram in eq. (A7) does not contribute at 
all in leading order, as can be seen by comparison with 
the first one, 



-—5V, h = — /i (^-^-^J = —h O (e) , 

(A9) 

where we have again used the power-counting rules. 
Thus, using eqs. (A7,A8,A9) gives 



T/(6) 

= ^ [II - '^hlt] (1 + 0(e)) + ^yi'- (AlO) 



The "nested" divergence proportional to /fimt^jdiv 
cancels as required and the two-loop counterterm and 
renormalized effective potential can be determined. They 
are given by the special case L = 2 of the general ansatz 



Veff{<t>)L 



y{2L+2) 

i-^) ^Ljr (A (^")) (1 + 0{e)) (All) 

V(2i+2) , 

___ (Jfinite(^„)^i^ (1 + Q^^^ (A12) 



Thus we obtain in leading order 

^effL = ^0 + ^ ^iV (1 + 0(e)) + Gl 



JV=1 



L-l „A2N+2) 

5Vl' + V L-t. /f(l + 0{e)) + 



N=l 



y(2L+2) 

2^L! 



5Vl' + 



y(2L+2) L 

2^L\ ^ I TV 

N=l ^ 



A (-A ) (l + ^^(e)) 



Y(2L+2) 

2^L! 



((/r*")^-(-/rT)(i+^^(e))- 

(A14) 

The last line shows that again all mixed terms in 
and /f'^ cancel, which is important for consistency, and 
implies the validity of the ansatz in eqs. (A11,A12) at 
order L, thereby completing the proof. Thus the final 
result for the renormalized second derivative of the ef- 
fective potential in bubble-approximation is given by the 
sum over all contributions from eq. (A12) {h= 1), 



E ^^^^^ (/r"=(V-"(0)))^ (1 + 0{e)). 

(A15) 



L=0 



Appendix B: HEAT KERNEL EXPANSION 



For L > 3, the upper eqs. can be proven by induction. 
Let us assume SV" {(j)) is given by eq. (All) for I < 
L — 1. Then we have to show eqs. (A11,A12) for L, 
which is done by evaluating the contribution of order 
to V^'ff as given in eq. (A3). This contribution can 
also be expressed as a sum of diagrams similarly to eqs. 
(A4,A7). Let G{N,l,ni,n2, . ■ . ,nL-i) denote the sum 
of all diagrams with A'' tadpoles, rifc insertions of SV^' 
into internal lines and with vertex y(2W+2) jf ; — g and 
Since the diagram should be of order 
h^, one has L = N+l + J^krik- This immediately implies 
that only counterterms of order less than L can enter 
for A^ > 1. All these contributions can be calculated 
using the ansatz (All). For this we note that the higher 
derivatives of 5V^' from eq. (All) have the form 



5Vt^ = (-l)'^I^ {ir{V")f (1 + 0(e)) (A13) 

in leading order in e, which can be seen in a similar way 

as in eq. (A8). Furthermore, using this relation, one 
can see that among the diagrams with fixed N and L all 
diagrams with counterterm insertions are suppressed by 
at least one factor of order e relative to Qj^ = G{N, L — 
N,0, . . . ,0), as in the two-loop case, see eqs. (A7,A9). 



The one-loop contribution to the effective action for 
the Lagrangian (40) is given by the determinant 

r[(/.]ii = -^lndeti, (Bl) 

with the operator A = n + X &ndX = d'^V{(t>, R)/d(fP'. 
The generalized zcta- function for A is Ca{v) = Ylm-^ni 
where are the eigenvalues of A. Using zeta-function 
regularization (see e.g. [62, 69]) the determinant can be 
written as 

r[<^]iL = |(C;i(0) + CA(0)lnM'), (B2) 

where we introduced a renormalization scale /x. The 
zeta-function can also be expressed via the heat ker- 
nel K(x,y,s) fulfilling the heat equation i-^K{x,y, s) = 
A{x)K{x, y, s) with K{x, y, 0) = ^(a; — y) as 

Ca{'^) = TT^ I ds{isY~^ I d^xK{x,x,s). (B3) 
J- Jo J 

The ansatz for K of Refs. [65, 70] is 
K{x,y,s) = j^;^^--^ A\fli{x,y)G{x,y,s) 

(B4) 
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where ij{x,y) is the proper arclength along the geodesic 
from X to y and Ay^ the Van Vleck-Morette de- 
terminant fulfiUing Avm{x,x) — -~g{x). Inserting 
this ansatz together with the expansion G{x,y,s) = 

Y^^.oii^VdA^^ y) i^to ^1- (^3) ^^'^ using eq. (B2) yields 
for the effective action 



IL 




where we have set X = X — R/Q. The coincidence limits 
of the gj can be calculated recursively. We quote the 
result for the lowest orders from Ref. [65], 

9o{x,x) = 1, gi{x,x) = 0, 

g2{x,x) = T^{R^,p„^R^''P^-R^,R^p-^DR+^DX 
= TM*^^ 360^^ 30'-'^+ e'-'^' 

(B6) 

where C and G are the Weyl- and Gauss-Bonnet terms 
as given in section IV. The gj with j > 3 contain higher- 



order curvature scalars built from the curvature- and 
Ricci tensors and space-time derivatives of R and X and 
correspond to finite contributions to the one-loop effec- 
tive action (B5), whereas the j = 0, 1, 2-contributions 
come along with divergences proportional to goX^, giX 
and g2- Using eq. (B6) one can see that it is neces- 
sary to introduce counterterms proportional to X^ = 
{d^V/d(l)^ - R/Qf, n{X - R/5) = a{d^V/d(f)^ - R/5), 
C and G in order to c;ancel these divergences, which is 
already done implicitly in the result (B5) for the effec- 
tive action through the zeta- function regularization [69] . 
Nevertheless, all opcirators contained in the counterterms 
should be already present in the tree level action. 
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